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Abstract 

Aims. The purpose of this study is to examine the prediction that the deuterated Hj ion, H2D+, can be found 
exclusively in the coldest regions of molecular cloud cores. This is also a feasibility study for the detection of the 
ground-state line of orriio-H2D+ at 372 GHz with APEX. 

Methods. The (lio — > In) transition of H2D^ at 372 GHz was searched towards selected positions in the massive 
star forming cloud OriB9, in the dark cloud L183, and in the low- to intermediate mass star-forming cloud R CrA. 

Results. The line was detected in cold, prestellar cores in the regions of OriB9 and L183, but only upper limits 
were obtained towards other locations which either have elevated temperatures or contain a newly born star. The H2D^ 
detection towards OriB9 is the first one in a massive star-forming region. The fractional ortho-HiD* abundances 
(relative to H2) are estimated to be ~ 1 10"'" in two cold cores in OriB9, and 3 10 '" in the cold core of L183. 

Conclusions. The H2D^ detection in OriB9 shows that also massive star forming regions contain very cold 
prestellar cores which probably have reached matured chemical composition characterized, e.g., by a high degree of 
deuterium fractionation. Besides as a tracer of the interior parts of prestellar cores, H2D+ may therefore be used to 
put contraints on the timescales related to massive star formation. 

1 Introduction 

The trihydrogen ion, H3 , is supposed to become the principal carrier of positive charge in the centres of cold, dense 
cores when 'heavy' elements like C, O and N, are nearly completely depleted (fW almsley et al. 2004| ). Because 
deuterium fractionation reactions are favoured in cold gas, relatively large abundances of the isotopoloques H2D^ 
and D2H^ are to be expected in these objects. This has been also confirmed by observations (ICaselli et al. 2003 1; 
IVastel et al. 20041 1. While infrared absorption spectroscopy of H3 can be used to extract vital information on the 
gas columns towards infrared sources (e.g. |McCall et al . 19991 ). the rotational lines of H2D^ and D2H^ probe the 
physical conditions of cold, prestellar cores. The abundance of H3 and its deuterated forms depend on the cosmic ray 
ionization rate of H2, and on the abundances of destructing agents: electrons, gaseous neutral species like CO and N2, 
and negatively charged grains (e.g. ICaselli et al. 2003 1; [Walmsle y et al. 2004| ). Furher more, the H3 abundance and 
the HjD' ^fli^ abundance ratio depend heavily on the ortho:para ratio of H2, which in turn is a function of time and 
density ( llPineau des Forets et al. 19911 : Flower et al 20061 1. 

The ground-state (1 10 — > 1 11) transition of ortho-H2D^ lies between adjacent atmospheric O2 and H2O absorption 
lines and its observation requires extremely good conditions. A reasonable limit is that the precipitable water vapour 
content of the atmosphere, PWV, is less than 0.5 mm, which can be achieved at high-altitude observatories only. H2D^ 
has been previously detected from Mauna Kea towards a protostellar core (ISt arketal , 19991), and in a small number 
of prestellar dark cloud cores (ICaselli et al. 2003|; jStarket al. 2004|; IVastel et al. 2006|). 

In this Letter we report on the first H2D^ observations with the Atacama Pathfinder Experiment, APEX, during 
its Science Verification periods in July, October and November 2005. The main goal of this series of observations 
was to test the feasibility of the 372 GHz line observations with this instrument. In the course of these measurements 
H2D^ was detected in a core belonging to a high-mass star forming region. This may open new vistas to the chemical 
evolution preceding the collapse of massive stars. 



'This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX). APEX is a collaboration between the Max-Planck- 
Institut fiir Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory. 



2 Source selection 



The selection of targets contains five starless cores of molecular clouds and one massive, cold core which encloses a 
low-luminosity far-infrared source. 

The massive, cold core OriB9 embedding IRAS 05405-01 17 in the region of Orion B is described in ICaselh & Myers 1994| , 
fCaselli & Myers 1995 1, and [Harju et al. 1993 1. The average ammonia linewidth is only 0.29 kms ' . The core does 
not stand out on the '''CO and C"^0 maps of ICaselh & Mye rs 1995 1, probably because of CO depletion. In this 
survey we have included the three N2H^(1 - 0) peaks found by |Caselli & Myers 1994|. These are likely to pinpoint 
separate clumps within the core. The clump associated with the IRAS source is likely to represent an early stage of 
collapse at which newly born stars have not yet disturbed their surroundings. The subsidiary clumps (OriB9 E and N) 
may be in a still earlier, pre-collapse phase. Depending on the time spent in the pre-collapse phase, chemical evolution 
may have resulted in a high degree of depletion and an increased H2D^ abundance. 

The two positions observed towards the nearby, starless dark cloud L183 (L134N) (see e.g. [iP agani et al. 20051 
and references therein) correspond to the 450/im continuum peaks in the SCUBA map of I Kirk e TaT'^OOSr ' The 
southern maximum (L183-S, the 850//m emission peak) can be assigned to a very cold, dense core with a high 
degree of molecular depletion. The two sources may represent different evolutionary stages of prestellar cores 
dLehtine n et al. 20031 V Therefore, it is interesting to compare their chemical and dynamical properties. After per- 
foming these observations we learned that the H2D^ has been detected at several positions along the dense ridge of 
LI 83 with the CSO by IVastel et al. 20061 . 

The starless, dense clump R CrA NW in the northwestern part of the R Coronae Australis cloud lies close to the 
1 .2 mm continuum source "MMS 10" (IChi ni et al.2003l ). Our position corresponds to an ammonia peak (Harju et al., 
in preparation). In this clump the kinetic temperature is higher than in the rest of our sample, probably implying that 
depletion is less marked. 
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' Tdust derived from the 450/850/im ratio (Kirk & Ward-Thompson 2006, private communication) 



3 Observations and data reduction 

The observations were made with APEX in July 2005 (RCrA and L183) and in October and November 2005 (OriB9). 
The 372421.364 MHz hne of H2D+ was observed in the upper side band with the APEX-2A SIS DSB receiver. The 
HPBW of the telescope is about 17" at this frequency. The observing mode was position switching with the off position 
-30' away in R.A. The integration time for each scan was 20 seconds. A calibration measurement was done every 
10 minutes. The first observations (RCrA) were made with the ASC 2048 channel autocorrelator using a bandwith 
of 128 MHz. Because the performance of the system was unsatisfactory with the ASC, the MPIfR Fourier transform 
spectrometer, EFTS, was used for all subsequent observations. The 1 GHz bandwith of the EETS was divided into 
16384 channels resulting in a channel width of 61 kHz which corresponds to ~ 49 ms ' at the observed frequency. 
The observing conditions ranged from excellent (PVW 0.2 mm, zenith opacity 0.24 at 372 GHz) to reasonable (PVW 
0.7 mm, zenith opacity 0.6 at 372 GHz). Depending on the weather and the elevation of the source, the DSB system 
temperature was between 130 and 300 K. 

Most of the observed spectra have ripple due to variations of sky emission, reflections in the telescope optics 
and instability of the receiver In the data reduction the possible low frequency ripple was first fit with a sinusoidal 
baseline whereafter possible higher frequency ripple was removed by masking the corresponding frequency in the 
Fourier transform. Finally, a first order baseline was subtracted around the source velocity. The mirror sideband of the 
receiver, centred at about 360 GHz, lies at a more transparent frequency than the signal sideband. The difference in the 
atmospheric opacity between the side bands was estimated using an atmospheric model and was taken into account in 
the calibration at the telescope. The telescope time spent on this project is 19.5 hours. 
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Figure 1: The observed H2D^(li() - In) spectra after Harming smoothing (Av = 98 ms Top: L183-S and L183-N. 
Middle: OriB9-E, OriB9-N and IRAS05405-0117. Bottom left: R CrA NW (ASC spectrometer). Bottom right: 
The probable detection of N2H+(4 - 3) at 372.67 GHz towards IRAS 05405-0117. The line can be seen in the same 
FFTS spectrum as the 372.42 GHz H2D+ line. 



4 Results 

The Hanning smoothed spectra are shown in Fig.^ A summary of the observations with the line parameters from 
Gaussian fits is presented in Table|2l 

The H2D+ Une is detected towards L183-S and 20" north and south of it. The line in L183-S (0,0) with the best 
S/N is single peaked, but suggests slight asymmetry. The FWHM (0.42 kms"') is a little larger than expected from 
thermal broadening at Tkin = 7 K (0.28 kms '). Our L183 positions lie near positions included in the N-S oriented 
H2D^ strip observed at the Caltech Submillimeter Observatory (CSO) by IVastel et al. 20061 . Their position Ad = 0" 
corresponds to the off'set (-4", +19") from L183-S. The antenna temperatures measured at APEX towards L183-S 
(0", 0") and (0", +20"), and the 3 a upper limit < 0.4 K obtained towads L183-N (close to Vastel's A6 = +80") 
are consistent with the CSO data. On the other hand, the line intensity at L183-S (0", -20"), T* ~ 0.9 + 0.2 K, is 
surprisingly high in view of the fact that Vastel et al. obtain antenna temperatures of x 0.7 K and 0.4 K towards 
the offsets (-4", -11") and (-4", -31") from L183-S, respectively. This suggests that the H2D^ distribution peaks 
slightly south of the dust emission peak but falls very steeply towards the southern end of the dense ridge. 

In OriB9, weak H2D^ lines with x 0.12 K at wlsr ~ 9.1 - 9.2 kms ' are detected towards the two starless 
condensations, OriB9 E and N. The LSR velocities agree with those of the previously observed NH3 (9.2 kms ' at 
OriB9 E and N) and N2H^ (9.2 kms ' at the IRAS position) lines. An upper limit of 0.10 K (3cr) is obtained towards 
IRAS 05405-01 17. This spectrum has, however, another line with ~ 0.4 K at an about 250 MHz higher frequency. 
The most probable identification is N2H+(4 - 3) at 372.67 GHz in the signal band (USB, ^lsr = 9.2 kms '). We 
estimate from the N2H+(1 - 0) data of [Caselh & Myers 1994| that the N2H+ column density towards IRAS 05405- 
0117 is ~ 6 - 8 10'^ cm"^. The intensity of the supposed N2H^(4 - 3) line is consistent with this value. This line is 
not detected in other spectra. 

Towards RCrA NW we obtained a 3cr upper limit of 0.3 K. The observations were made at low elevations. Further- 
more, the system temperature was higher than expected from the observing conditions, probably because of problems 
with the integration of the ASC correlator into the system. 



5 Column densities and abundances 



The observed H2D^ line has quadrupole hyperfine structure due to the spins of D and the two H nuclei. The splitting 
is, however, very small (Av ~ 80 kHz, IJensenetal. 199l1 ) compared with the Doppler width at 10 K (Av ~ 420 kHz), 



Table 2: Line parameters derived from Harming smoothed spectra with a velocity resolution of 98 ms 
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and we treat the line as if it had a single component. We estimate the ortiio-H2D'^ column densities, N{o — HaD^), in 
the same manner as done in ICaselli et al. 2 003 1. A lower limit for the excitation temperature, Tex, is obtained from 
the observed by assuming that the line is optically thick. The kinetic temperature of the gas, Tkin, sets an upper 
limit. 

Towards L183-S (0,0), the reasonable Tex range is 6.3 K < Tex < lOK (allowing Tkin be slightly higher than the 
dust temperature). Assuming that T^.^ = T^^^st = 7 K, and that the source fills the main beam uniformly, we arrive at the 
values To = 1.25 and A^(o-H2D+) = 2.7 10'^ cm~^. The main beam efficiency, jjmb, is assumed to be 0.7. This column 
density is very close to that found by DCaselli et al. 2003ll towards the centre of LI 544. The peak H2 column density 
towards L183-S derived from SCUBA data is 9.1 10^^ cm"^ (Kirk & Ward-Thompson, 2006, private communication). 
Using this value we obtain the fractional ortho-H2D^ abundance X{o - H2D^) = 3 10 By varying Tex in the given 
range, N{o - H2D^) and X{o - H2D^) decrease (higher Tex, smaller tq) or increase (lower Tex, larger tq) by a factor of 
three. The line shape does not support, however, the idea of large opacities, and we think the values corresponding to 
Tex = 7 K are the most likely. 

OriB9 is little studied and dust continuum measurements are not available. The Tkin derived from ammonia is 
10 K, and the minimum Tgx from the line intensity is 4 K. Using the assumption Tex = 7 K, which is midway between 
the two extremes, we get tq ~ 0.13 and N{o - H2D^) ~ 3.0 10'^ cm"^ towards N and E. The ammonia column density 
in these positions is ~ 10'^ cm"^. Assuming that the fractional NH3 abundance is 3 10"^ (e.g. |Harju et al. 1993| ) 
we get X{o - H2D^) ~ 1 10 The line profiles with rather poor S/N ratios do not exclude large optical thicknesses. 
The column density range implied by the possible Tex values is ~ 1.0 10'^ - 5 10'^ cm"^, and the derived fractional 
abundance has the corresponding uncertainty. 

According to the model of |Walmsley et al. 2004| , the characteristic steady-state value of the o/p-ratio of H2D^ is 
~ 0.3 in the density range ~ 10^ - lO*" cm"-* appropriate for the objects of this study. Adopting this o/p-ratio the 
total H2D+ abundances in L183-S and OriB9 become ~ 1.3 10 and ~ 4 10 

6 Discussion 

The performance of APEX and its equipment, and the atmospheric transmission on Chajnantor are found to meet 
very high standards. The telescope is therefore likely to become a very important tool for studies of molecular cloud 
interiors and star formation using H2D^, and other 'difficult' molecules. 

The present observations towards a small sample of dense cores with some diversity of physical characteristics 
suggest that either an elevated temperature (as in R CrA NW) or the presence of an embedded star (as towards IRAS 
05405-01 17), even if there is little evidence for star-cloud interaction, decreases the changes to find H2D^. The kinetic 
temperature, velocity dispersion and the fractional H2D^ abundance in OriB9 are similar to those in the prestellar dark 
cloud cores L1544, 16293E, and LI 83, where strong emission of this line has been detected previously. The masses and 
the central densities of these nearby cores are 2-3 Mq and ~ 10^ cm"^, respectively ([ Vastel et al. 2006 1 and references 
therein). The total mass of the OriB9 core estimated from ammonia is of the order of 100 Mq ( |Harju et al. 1993) ), 
and the subcondensations seen in the N2H^ map (fCaselli & Myers 1994]) are likely to be an order of magnitude more 
massive than L1544, 16293E, or L183. A (sub)millimetre continuum map is needed to confirm this. Nevertheless, 
OriB9 seems to be capable of forming a massive star or a dense cluster of low- to intermediate-mass stars, and is 
therefore exceptional among sources detected in H2D^ so far 

This detection confirms the existence of very cold, quiescent, dense cores in massive star forming regions. The 
previous evidence for such objects is scarce. Some infrared dark clouds (IRDCs) have gas temperatures approaching 



10 K and ammonia linewidths slightly below 1 kms"' ( IPillai et al. 20061). Yet another massive core with these 
characteristics has been recently found in the region of ISOSS J18364-0221 (|Birkmann et al. 2006 1). The linewidths 
in these objects are, however, clearly larger than in OriB9. On the other hand, only a small fraction of IRDCs have 
been observed in spectral lines to date. Because compression leads to an intensified cooling by molecules and dust, the 
collapse of all dense cores should be preceded by a cold, quiescent phase. Very cold GMC cores may have remained 
indiscernible because of their short life-time or the fact that large-scale surveys are usually biased towards the presence 
of a certain molecular species which might be depleted in the coldest regions. 

As discussed recently by | Flower et al 2006 1, the HiD^ abundance depends inversely on the ortho:para ratio of 
H2, which is largest at early stages of core evolution. Consequently, a high degree of deuterium fractionation is a 
sign of matured chemistry characterized by a low o/p H2 and probably a high degree of molecular depletion. The 
ortho-H2D^ detection towards OriB9 suggests an evolved chemical stage and tells of a longlasting prestellar phase. 
However, estimates of the o/p ratio of H2 and the degree of deuterium fractionation using other tracers are needed to 
confirm this. Despite these obstacles, it seems viable to use the 372 GHz H2D^ line together with chemistry models 
to estimate timescales related to the early evolution of massive cores. 
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